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Abstract. - Recent advances in nanofabrication and optical control have garnered tremendous 
interest in mult i-qubit- cavity systems. Here we analyze a spin-glass version of such a nanostruc- 
ture, solving analytically for the phase diagrams in both the matter and radiation subsystems in 
the replica symmetric regime. Interestingly, the resulting phase transitions turn out to be tunable 
simply by varying the matter-radiation coupling strength. 



Atomic physics, nanostructurc materials science and 
quantum optics have recently made remarkable advances 
in the fabrication and manipulation of matter-radiation 
systems [1]. The energy gap in a semiconductor quan- 
tum dot [2] can be engineered by varying the dot size and 
choice of materials. For example, vanishingly small optical 
gaps could be obtained using InAs/GaSb or HgTe/CdTe 
quantum dots, while vanishingly small inter-subband gaps 
can be obtained by increasing the dot's size [2]. Hence 
tailor-made two-level 'qubit' (i.e. quantum bit) systems [3] 
can be built on the length-scale of 10^ — 10"^ Angstroms, 
and with various geometric shapes, using a range of III-V 
and II- VI semiconductors [2]. In addition, experimental 
control of the qubit-cavity coupling A has already been 
demonstrated [4-6] for quantum dots coupled to photonic 
band-gap defect modes [7], as well as for atomic and su- 
perconducting qubit systems. A qubit-qubit interaction 
can arise, for example, from the electrostatic inter-dot 
dipole-dipole interaction between excitons and/or conduc- 
tion electrons, and can be engineered by adjusting the 
quantum dots' size, shape, separation, orientation and 
the background electrostatic screening. The interaction's 
anisotropy can be engineered by choosing asymmetric dot 
shapes. Disorder in the qubit-qubit interactions can be 
introduced by varying the individual dot positions dur- 
ing fabrication, or will arise naturally for self-assembled 
dots [2,8]. All the pieces are therefore in place for engi- 
neering all-optical realizations of condensed matter spin- 
based systems. 

Given these exciting developments in multi-qubit-cavity 
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nanostructurc, we study here the effect of a photon field 
on a set of qubits with disordered interaction, which can 
be viewed as a spin-glass [9, 10]. More specifically, we 
provide an analytic analysis of the phase behaviour with 
respect to temperature, spin-spin coupling strength and 
photon-spin coupling strength variations. We find that 
phase transition phenomena arise within both the spin 
(i.e. matter) and boson (i.e. photon) subsystems. With 
the current technology, an order of 10"^ quantum dots can 
be embedded in a cavity structure [5], the phase transi- 
tion should therefore be a prominent feature of the system. 
Also importantly from an experimental point of view, the 
resulting phase diagrams can be explored within a given 
nanostructurc array by varying the qubit-cavity coupling 
strength A. Furthermore, single quantum dot readout is 
currently under intense development (e.g., see [11]). This 
opens up the possibility of studying the local feature of 
the system, which would be extremely helpful in under- 
standing spin glasses. 

Our main results follow from solving analytically, in 
the replica symmetric regime, a generalised version of the 
Dicke model [12,13]. The Dicke model was originally de- 
veloped to describe the radiative decay of a gas of two- level 
systems [12], the superradiance-subradiance phase transi- 
tion was later discovered [14]. Debates then ensued as 
to whether the phase transition is physical, due to the 
constraint of the Thomas-Reiche-Kuhn (TRK) sum rule 
(e.g., see [15]). However Keeling has recently provided 
a convincing argument for the physical existence of the 
Dicke phase transition by analysing the full atom-photon 
hamiltonian [16]. Furthermore, in contrast to the debates 
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concerning atom-photon interactions, our primary concern 
here is solid-state systems such as quantum dots where 
there are usually many electrons. As such, the dipole 
strengths can be re-distributed among the different en- 
ergy levels, and hence the constraint from the TRK sum 
rule can be drastically modified for the lowest two energy 
levels ^. For these reasons, we believe that it is indeed 
legitimate for us to employ here the original Dicke hamil- 
tonian, supplemented with a spin-spin coupling term to 
model our qubit-qubit interaction. More specifically, our 
hamiltonian is of the form: 

^ A e 
^ = + E + + ^7) + E l^f + 

i=i j 

(1) 

where the operators a, and , correspond to the 
photon field and quantum dot j respectively, and Hss is 
the spin-spin coupling term. 

We now make the following assumptions: 

1. Hss is of the form J2i<j Jij<^f<^f ■ 

2. e< min{l,A///V,{Jy}} 

As discussed in [6] , there arc many ways to realize the first 
assumption. For example, each quantum dot can be engi- 
neered to have an elongated form along the a;-direction, by 
biasing the growth process along this direction. Applying 
an electric field along x, will then create large permanent 
dipole moments in that direction. One can use undoped 
dots, in which case the dipole results from the exciton, 
or doped dots, in which case the dipole originates from 
the conduction-subband electron biased along x. The first 
assumption also implies that the coupling in the x direc- 
tion overwhelms the (multipole) coupling in the y and z 
directions. The second assumption requires a negligible 
energy gap between the two levels of the quantum dot. 
This could be achieved using HgTe/CdTe quantum dots 
in order to reduce e, and/or by engineering strong dot-dot 
interactions { Jy} and a strong cavity-dot optical coupling 
A. 

We would like to note that we have studied photon-spin- 
glasses systems in a more general setting in [17] where the 
phase diagram is analysed by numerically solving a set of 
self-consistent equations deduced using the Trotter-Suzuki 
method [18]. Although the hamiltonian employed here is 
more restrictive by comparison, it has the virtue of being 
amenable for analytical treatment as we shall see below. 

For self-assembled dots, the dot-dot coupling terms Jij 
will have an inherent disorder. Alternatively, such disor- 
der can be built in during growth by varying the dot-dot 
separations. As a general model for disorderness, we make 



the assumption that the distribution of Jij 's is Gaussian: 



^Interestingly, the borrowing of dipole strength from other energy 
levels seem to be exploited by some photosynthetic systems in nature 
as discussed by Hu et al, Quart. Rev. Biophys. 35 (2002) 1. 



(2) 



with Jo and J representing the mean and standard devi- 
ation of the probability distribution. We note that nega- 
tive Jij is experimentally feasible given the multipole (e.g. 
dipole-dipole) nature of the qubit-qubit interactions. 

We now introduce the Glauber coherent states |a) [19], 
which have the following properties: a la) = a\a), {a\a^ = 
{a\a*, and / |q:)(q;| = l. In terms of this basis, 

the canonical partition function can be written as: 



Z{N,T) 



dRje{a)dlm{a) 



(ajTre" 



-fJH 



\a) . (3) 



We adopt the same assumptions as in [14]: (i) the order 
of the double limit in limAr^oo limfl_»oo X]^=o ^h"'* '^^^ 
be interchanged, and (ii) a/^fN and a) /^/N exist as iV ^ 
oo. With these assumptions, we rewriting Z{N, T) in the 
following form: 



Z{N,T) = j ^^e-^l«l'TVe-'5^' 



(4) 



where 



^' = E^-f-E^^.-f-f • (5) 

3 i<3 

It should now be clear why we made the assumptions con- 
cerning the form of our hamiltonian all the terms in H' 
are commutative and so we can integrate out a analyti- 
cally. Again, if the terms are not commutative, numerical 
methods will have to be used [17]. 

Performing the Gaussian integral with respect to a, we 
obtain 



Z = -^Trexp 



i<j 



2A^ 
N 



+ /3A2 



(6) 

Hence the problem can be mapped onto the traditional 
spin-glass hamiltonian if we make the transformation 
Jo Jo in the Jij probability distribution, where Jo = 
(Jo -I- 2A^). In other words, the phase diagram for the 
matter system (e.g., a nanostructure array) is equivalent 
to the usual spin-glass one [10] (c.f. Fig. 1). 

We now turn to consider the supcrradiant and subra- 
diant photon states in the optical subsystem. We recall 
that the order parameter for the subradiant-superradiant 
phase separation is defined to be [14] : 



e = 



(7) 



In view of this, we insert a factor h in the exponent — /3|a|^ 
in Z, i.e.. 



(8) 



p-2 



Spin-glasses in optical cavity 



Fig. 1: Phase diagram of the matter (e.g. quantum dot) subsys- 
tem within a multi-qubit-cavity such as a disordered quantum 
dot array coupled to an optical cavity mode. The scaled pa- 
rameter Jo can be varied simply by changing the qubit-cavity 
coupling strength A (see text). T is the temperature. 
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Fig. 2: Phase diagram showing the superradiant-subradiant 
transition in the optical subsystem. In the absence of spin- 
spin interactions, the system is always superradiant due to the 
two-level system's negligible energy gap e. 
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so that 



1 d[F] 



N dh 



(9) 



h=l 



We now integrate out a in Z and we obtain 

2A2 



hN 



(10) 

Namely, the only modification is the extra h in the de- 
nominator in the underlined term. 

We now recall the definition of the free energy for spin 
glass systems [10,20]: 

[F] = -i[logZ] = d[J]P(J,,) log Z (11) 

where d[J] = Ylij dJij and we employ the replica method 

to approximate [logZ] as lim„^o — ^ — ■ Going through 
the standard procedure of integrating out the quenched 
disorder in Js [20], we obtain: 



"1 — „Nf3^J^n/4 



[Z"] = 



d[q]d[r 



(12) 



where d[q] = n„<^ dquv, d[m] = n„ dniu and 
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uv 
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with 



(13) 



(14) 



where u, v are the replica indices and the X-superscripts 
in the as are dropped for clarity. Note that Jq is now 



(Jo + 2\'^ /h) and so we obtain for d[Z"]/dh the following 
expression: 



d[g]d[m]|^^^^(-m„ + 2(a")L)m„|e" (15) 



where (ct")l = ^"^^l ■ It follows that 

1 d[F] 1 d[Z"] 



N dh 
2A2 



h=l 



n[3N dh 



h=l 



j d[q]d[m] |E(2('^")i - mu)m}^ . (16) 

Since is proportional to N ^ we can evaluate the integral 
by steepest-descent. In the thermodynamic limit N ^ oo, 
we find that 



9 = const. X — y^(2((T") L — "^n)"^n^ 



(17) 



Since 51 is optimized when 
1 d 



qu 



and 



/3V2 ^qu^ 

J d_ 

13 Jq dniu 



logTre^ = (aV")L (18) 



logTre^ 



(19) 



we have 9 oc — m^. Therefore, in the replica- 
symmetric case, 9 cx and one can then derive the 
superradiant-subradiant phase diagram as shown in Fig. 
2, outlining the region where m ^ 0. Finally we note, by 
observation of Eq. (7), that spin-glass behaviour can also 
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arise in a multi-qubit-cavity system with disorder in both 
{A} and {J}, as long as some of the J's are negative. 

In conclusion, we have analyzed a novel optical realiza- 
tion of a Hamiltonian system which is of great interest 
within the condensed matter community, and have de- 
duced analytically the corresponding phase diagrams. In 
contrast to traditional realisations using magnetic solids, 
the phase transitions in this system can be explored simply 
by changing the matter-radiation coupling strength. 
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